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METHOD AND DEVICE USING MYOELECTRICAL ACTIVITY FOR 
OPTIMIZING A PATIENT'S VENTILATORY ASSIST 



5 FIELD OF THE INVENTION 

The present invention relates to a method and device for determining a 
level of ventilatory assist to a ventilator-dependent patient. 

BACKGROUND OF THE INVENTION 

10 Both the tension developed by a patient's muscle [34] and the duration of 

the muscle contraction [2] are factors that lead to respiratory muscle fatigue; 
these two factors can be expressed by indices such as the tension-time index [3] 
and the pressure-time product [10, 20, 32, 35]. Bellemare and Grassino [3] 
showed a direct inverse relationship exists between the time of endurance of a 

1 5 fatiguing diaphragm contraction and the rate of decay of the ratio of the high to 
low spectral components (H/L) of the electrical activity EAdi of the patient's 
diaphragm, indicating that these two values are indicative of progressive failure 
to sustain load. The force exerted by the muscle has been shown to be directly 
related to the rate of decay of the power spectrum center frequency or the rate 

20 of decay of the above mentioned ratio H/L, and the level at which this power 
spectrum center frequency or ratio H/L plateaus [16, 21, 28]. Such shifts in the 
power spectrum reflect a reduction in the muscle action potential conduction 
velocity [28, 38, 39], and constitute an early indication that, at the cellular level, 
these breathing patterns cannot be maintained indefinitely [3]. 

25 

Hyperinflation, which impairs the length-tension relationship of the 
respiratory muscles, i.e. the transformation of the neural activation into a 
mechanical output or pressure, reduces the capacity of the respiratory muscles 
to generate pressure (neuromechanical uncoupling), unless the electrical activity 
30 EAdi of the patient's diaphragm is increased. Studies have shown that when the 
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inspiratory pressure, flow and duty, cycle remain constant, increases in end- 
expiratory lung volume (EELV) promote reductions in endurance time [33, 44] 
and sustainable pressure [11]. In an animal model, Tzelepis et al [44] proposed 
that, under these conditions, diaphragm shortening would require greater 
5 excitation to generate a given sub-maximum tension, and that this increased 
excitation might account for the greater contractile muscle fatigability observed 
at shorter muscle length. 

The level of partial ventilatory assist, with the aim to ensure adequate 
10 pulmonary ventilation while preserving inspiratory muscle function, is generally 
set on an empirical basis in the clinical setting. 

It has been proposed that an optimal level of partial ventilatory assist 
could be determined from the lowest stable breathing frequency f B achieved, i.e. 
15 without bradypnea or apnea. In patients, this corresponded to 16.4 bpm (breaths 
per minute) and was associated with a tidal volume V T of 11.8 ml/kg. However, 
mechanical lung modeling in that study demonstrated that such a level of 
support actually resulted in a total unloading of the respiratory muscles. 

20 Others have defined an optimal level of partial ventilatory assist as that 

which produces the lowest swings of transdiaphragmatic pressure P d j and found 
that this condition was associated with a breathing frequency f B of 19.7 bpm and 
a tidal volume V T of 11 .7 ml/kg. The transdiaphragmatic pressure P d i in the latter 
study was used as a marker of inspiratory effort. 

25 

Jubran et al [20] defined an upper bound inspiratory pressure-time 
product lower than 125 cm H 2 Os/min as a desirable level of inspiratory effort to 
be achieved during partial ventilatory assist. Although arbitrarily determined, this 
threshold was justified by the fact that it corresponded to a tension-time index 
30 TTdi well below that considered to indicate impeding inspiratory muscle fatigue. 
The study found a high variability in pressure-time products between patients 
and demonstrated that a breathing frequency f B < 30 bpm and a tidal volume V T 



WO 2005/077268 PCT/CA2005/000217 



of 0.6 L were better determinants of an optimal level of inspiratory effort during 
partial ventilatory assist. Based on these breathing pattern findings, it is likely 
that the level of respiratory muscle unloading provided by this method of 
optimizing partial ventilatory assist was lower than that of the above discussed 
5 studies. 

Brochard et al [8] defined an optimal partial ventilatory assist level as the 
lowest level of ventilatory assist, which when implemented, maintained the 
highest level of diaphragmatic electrical activation without the occurrence of 
10 fatigue as evaluated via power spectrum analysis of the electrical activity EA d i of 
the patient's diaphragm. Interestingly, such levels of partial ventilatory assist 
were associated with a breathing frequency f B of 20-27 bpm and a tidal volume 
V T of 8.0 ml/kg, these values being similar to those later reported by Jubran et al 
[20]. 

15 

SUMMARY OF THE INVENTION 

In accordance with the present invention, there is provided a method for 
20 determining a level of ventilatory assist to a ventilator-dependent patient, 
comprising: calculating a critical threshold of a respiration-related feature, 
wherein fatigue of a respiratory muscle of the ventilator-dependent patient 
develops when the critical threshold is reached by the respiration-related 
feature; and controlling the level of ventilatory assist to the ventilator-dependent 
25 patient in relation to the critical threshold of the respiration-related feature so as 
to prevent fatigue of the patient's respiratory muscle. 

The present invention also relates to a device for determining a level of 
ventilatory assist to a ventilator-dependent patient, comprising: a calculator of a 
30 critical threshold of a respiration-related feature, wherein fatigue of a respiratory 
muscle of the ventilator-dependent patient develops when the critical threshold 
is reached by the respiration-related feature; and a controller of the level of 
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ventilatory assist to the ventilator-dependent patient in relation to the critical 
threshold of the respiration-related feature so as to prevent fatigue of the 
patient's respiratory muscle. 

The foregoing and other objects, advantages and features of the present 
invention will become more apparent upon reading of the following non- 
restrictive description of illustrative embodiments thereof, given by way of 
example only with reference to the accompanying drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 
In the appended drawings: 



Figure 1 is a schematic representation of a non-limitative example of 
15 experimental set-up for measuring diaphragm's electrical activity EAdi, 
esophageal pressure P es , gastric pressure P ga , respiratory airflow and tidal 
volume Vt, and for displaying on line the target transdiaphragmatic pressure P d i 
and the root-mean-square (RMS) of the diaphragm's electrical activity EAdi; 

20 Figure 2 are illustrative examples of tracings of tidal volume V T , 

diaphragm electrical activity EAdi, transdiaphragmatic pressure P d j, esophageal 
pressure P es , and gastric pressure (P ga ) measured on a subject during "volume" 
maneuvers and "expulsive" maneuvers; 

25 Figure 3 are examples of bar graphs displaying drops in center frequency 

CFdi, targeted levels of transdiaphragmatic pressure Pdj, diaphragm pressure- 
time product PTPdi and the associated diaphragm's electrical activity EAdi 
observed during volume and lower-pressure expulsive and higher pressure 
expulsive maneuvers; 

30 

Figure 4 are examples of graphs from one representative subject 
showing the center frequency CF d j, the root-means-square (RMS) of the 
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diaphragm's electrical activity EAdi and the diaphragm pressure-time product 
PTPdi plotted over time during the volume maneuver (circles) and the two 
expulsive maneuvers at end-expiratory lung volume (EELV), one targeting a 
lower transdiaphragmatic pressure P di (squares) and the other a higher P d i 
5 transdiaphragmatic pressure Pdi (triangles); 

Figure 5 is a flow chart and block diagram of a first non-restrictive 
illustrative embodiment of the method and device according to the present 
invention, for determining a level of ventilatory assist to a ventilator-dependent 
1 0 patient; and 

Figure 6 is a flow chart and block diagram of a second non-restrictive 
illustrative embodiment of the method and device according to the present 
invention, for determining a level of ventilatory assist to a ventilator-dependent 
1 5 patient. 



DETAILED DESCRIPTION OF THE ILLUSTRATIVE EMBODIMENTS 

A study was conducted to determine in humans whether an increased 
20 electrical activity EAdi of a patient's diaphragm, with neuromechanical 
uncoupling, promotes greater reductions in the center frequency CF d i of the 
diaphragm's electrical activity EAdi, when the diaphragm pressure-time product 
PTPdi is kept constant. An additional aim of the study was to establish the extent 
to which the diaphragm pressure-time product PTP d j needs to be increased, in 
25 the presence of normal neuromechanical coupling, in order to reproduce the 
drop in center frequency CF d j observed with uncoupling. 

More specifically, the study evaluated whether increased diaphragm 
activation induced by an increased lung volume promotes increased drops in the 
30 center frequency CF d j of the diaphragm's electrical activity EAdi when the 
diaphragm pressure-time product PTP di is kept constant. Five healthy subjects 
performed runs of intermittent quasi-static diaphragmatic contractions with a 
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fixed breathing pattern. In separate runs, the subjects targeted 
transdiaphragmatic pressures P d j by performing end-inspiratory holds at total 
lung capacity with the glottis open (neuromechanical uncoupling), and at end- 
expiratory lung volume by performing expulsive maneuvers (no 
5 neuromechanical uncoupling). Diaphragm activation and pressures were 
measured with an electrode array and with balloons, respectively, mounted on 
an esophageal catheter. Reproduction of a transdiaphragmatic pressure P dj of 
«31 cm H 2 0 during neuromechanical uncoupling increased lung volume to 
77.5% of the inspiratory capacity, increased the diaphragm's electrical activity 

10 EAdi from 25% to 61% of the maximum and resulted in a 17% greater drop in 
center frequency CF d i. In order to reproduce, in the absence of neuromechanical 
uncoupling, the decrease in center frequency CF di observed during 
neuromechanical uncoupling, a two-fold increase in transdiaphragmatic pressure 
P d i and diaphragm pressure-time product PTP d j was required. It was concluded 

15 that a constant diaphragm pressure-time product PTP di does not necessarily 
result in a center frequency CF di of the diaphragm's electrical activity EAdi that 
remains stable when activation is increased. 

METHODS 

20 

Subjects 

Five healthy subjects (1 female, 4 males) with a mean age of 40.6 ± 8.0 
years participated in the study. The study was approved by the Scientific and 
25 Ethical Committees of Sainte-Justine's Hospital and all subjects gave their 
informed consent. 

Experimental protocol 

30 Figure 1 is a schematic representation of a non-limitative example of 

experimental set-up. On the left, an esophageal catheter-mounted multi- 
electrode array 5 is used to measure diaphragm's electrical activity EAdi and 
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balloons 8 and 9 mounted on the catheter on opposite sides of the electrode 
array 5 are used to measure esophageal pressure P es and gastric pressure P ga . 
The catheter 6 was passed trans-nasally and positioned at the gastro- 
esophageal junction 10. Respiratory airflow was measured with a 
5 pneumotachograph 3 and tidal volume V T was obtained by integrating inspiratory 
flow. On the right, the target transdiaphragmatic pressure P d i and the root-mean- 
square (RMS) of the diaphragm's electrical activity EAdi are displayed on line. 

Referring to Figure 1, each human subject 2, while seated in an upright 
10 chair (not shown) and facing a the monitor 1 of a computer 4, performed 
repeated maximal inspirations to total lung capacity (TLC) in order to obtain 
three reproducible voluntary maximum values for the diaphragm's electrical 
activity EAdi. Each subject 2 was subsequently asked to perform intermittent, 
near-isometric diaphragmatic contractions of 10 seconds duration, separated by 
15 5 seconds relaxation periods during which free breathing was allowed. With 
visual feedback of the transdiaphragmatic pressure P d i on the monitor 1 of the 
computer 4 a low level of transdiaphragmatic pressure P d i was targeted during 
two runs, while a higher level of transdiaphragmatic pressure P d i was targeted 
during a third run. The duty cycle was imposed by a sound signal, and each run 
20 lasted until a plateau in center frequency CF dj was reached, or until the subject 
was no longer able to maintain the target transdiaphragmatic pressure P di . 

Figure 2 are examples of tracings of tidal volume V T , diaphragm electrical 
activity EAdi, transdiaphragmatic pressure P di , esophageal pressure P es , and 

25 gastric pressure P ga measured in one subject during "volume" maneuvers and 
"expulsive" maneuvers performed during the hereinafter reported study. The 
"volume" maneuver consisted of an end-inspiratory hold at an increased lung 
volume, which resulted in the generation of a low P di (left tracing), whereas the 
two expulsive maneuvers were performed at end-expiratory lung volume 

30 targeting a lower P dj (middle tracing) and higher P di (right tracing). 



In order to obtain two different levels of diaphragm's electrical activity 
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EAdi for the same target transdiaphragmatic pressure P di , each subject 2 was 
. instructed to perform two different maneuvers: 

1. Volume maneuver: the subjects inspired close to their total lung capacity 
5 (TLC) and produced a given level of transdiaphragmatic pressure P d i (Figure 

2; left tracing). The transdiaphragmatic pressure P d i was maintained at this 
lung volume with the glottis open. 

2. Expulsive maneuver: the subjects performed expulsive maneuvers in order 
10 to generate a target transdiaphragmatic pressure P d j. All expulsive 

maneuvers were performed at end-expiratory lung volume (EELV) at lower 
and higher transdiaphragmatic pressures P dj (Figure 2, middle and right 
tracings). 

After having initially performed a volume maneuver run, each subject 2 
then performed two expulsive maneuver runs. One expulsive maneuver run 
targeted a transdiaphragmatic pressure P dj (lower pressure) similar to that 
observed during the volume maneuver but requiring less diaphragm's electrical 
activity EAdi, while another expulsive maneuver run targeted an increased 
transdiaphragmatic pressure P di (higher pressure) to reproduce the center 
frequency CF dj observed during the volume maneuver run. The volume 
maneuver was subsequently repeated once for retest purpose. The subject 
rested for 20 minutes between subsequent runs. 

25 Instrumentation 

Using the set-up of Figure 1 : 

- airflow and tidal volume were measured by a computer 4 through a 
30 pneumotachograph 3; 

- electrical activity EAdi of the patient's diaphragm was measured by the 



15 



20 
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computer 4 through the linear array 5 of electrodes mounted on an 
esophageal catheter 6 inserted through the patient's nostril (or patient's 
mouth) until the electrode array 5 is positioned in the gastro-esophageal 
junction 10 of the patient's diaphragm 7; 

5 

- esophageal P es and gastric P ga pressures were measured by the computer 4 
through the gastric 8 and esophageal 9 balloons mounted on the catheter 6 
on opposite sides of the array 5 of electrodes; and 

10 

- the transdiaphragmatic pressure P d j was obtained by the computer 4 by 
subtracting the measured esophageal pressure P es from the measured 
gastric pressure P ga . 

1 5 On-line automatic processing of diaphragm's electrical activity EAdi 

The diaphragm's electrical activity EAdi, more specifically a root-mean- 
square (RMS) EAdi signal was acquired, processed and displayed on-line using 
a standardized methodology [4, 36, 41]. The center frequency CF d i was 

20 evaluated for signal quality using established indices and criteria in accordance 
with a method disclosed by Sinderby et al [40]. To avoid influence of power 
spectral shifts on the EAdi signal strength, the RMS EAdi signal was calculated 
on the spectral moment of order 1 (M1) which is insensitive to conduction 
velocity [6] (see upper trace on the computer monitor 1. For more extensive 

25 review reference is made to Aldrich et al [1]. 

Off-line signal analysis 

Inspiratory duration Tj, total breath duration T to t, and breathing frequency 
30 f B , diaphragm's electrical activity EAdi and pressures P es and P ga were 
determined using the transdiaphragmatic pressure Pdi- The diaphragm pressure- 
time product PTP di was obtained by multiplying (i) the under-the-curve area 
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subtended by the P di signal by (ii) the breathing frequency f B . The amplitude of 
the signal of the diaphragm's electrical activity EAdi was expressed as a 
percentage of the voluntary maximum diaphragm's electrical activity EAdi 
obtained from TLC maneuvers [37]. Variables were compared between each of 
5 the maneuvers performed using one-way repeated measurements analysis of 
variance (ANOVA) and post hoc contrasts of significant effects were performed 
using the Student-Newman-Keuls test. Test-retest reliability of the P d1l EAdi and 
CFdi values obtained during the volume and expulsive maneuvers was evaluated 
by calculating the interclass correlation coefficient (ICC). 

10 

RESULTS 



The subjects were able to perform all maneuvers and maintain the 
imposed duty cycle (P=(X93; one -way ANOVA) during all protocols (Table 1). 

15 

Table 1 

Breathing pattern and targeted P di values during the three maneuvers performed 





Expulsive maneuver 


Expulsive maneuver 


Volume maneuver 








Lower pressure 


Higher pressure 


Subject 


Ti/Ttot 


Vt%IC 


Pdi 


Ti/Ttot 


Pdi 


Ti/Ttot 


Pdi 


1 


0.66 


89.2 


47.8 


0.66 


45.1 


0.67 


77.3 


2 


0.66 


65.6 


25.3 


0.67 


26.6 


0.66 


93.3 


3 


0.67 


68.3 


10.2 


0.65 


11.6 


0.66 


39.3 


4 


0.67 


88.8 


38.6 


0.66 


39.8 


0.64 


77.3 


5 


0.65 


75.5 


34.1 


0.67 


35.0 


0.67 


50.3 


Mean 


0.66 


77.5 


31.2 


0.66 


31.6 


0.66 


67.5 
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(±SD) 


(0.01) 


(11.1) 


(14.2) 


(0.01) 


(13.1) 


(0.01) 


(22.0) 



Values are means for each subject of all the maneuvers performed. Ti/Ttot, duty cycle; P d i, 
transdiaphragmatic pressure; Vt, tidal volume; IC, inspiratory capacity. All subjects were 
able to maintain the imposed duty cycle. 



Figure 3 are examples of bar graphs displaying drops in center frequency 
5 CFdj, targeted levels of transdiaphragmatic pressure P d j, diaphragm pressure- 
time product PTPdi and the associated diaphragm's electrical activity EAdi 
observed during the three volume and lower-pressure expulsive and higher 
pressure expulsive maneuvers performed in the study. The bars of the graphs of 
Figure 3 are average values obtained for the five subjects (± SD (Standard 
1 0 Deviation)). 

As shown in Table 1 and Figure 3, subjects were able to achieve and 
maintain similar target levels of transdiaphragmatic pressure P d i during the 
volume maneuver (high lung volume) and the lower-pressure expulsive 

15 maneuver at EELV. During the volume maneuver, subjects inspired to an 
average of 77.5±11.1% of their inspiratory capacity. In order to generate the 
same diaphragm pressure-time product PTP d j at different lung volumes, the 
volume maneuver (neuromechanical uncoupling) required a diaphragm's 
electrical activity EAdi of 60±8% of maximum compared to 25±8% for the 

20 expulsive lower-pressure maneuver at EELV. As shown in Table 2 and Figure 3, 
despite a matching of the diaphragm pressure-time product PTP d j, the volume 
maneuver promoted a 17% larger drop in the center frequency CF di than the 
expulsive low-pressure maneuver at EELV. Figure 4 are examples of graphs 
from one representative subject showing the center frequency CF d j, the root- 

25 means-square (RMS) of the diaphragm's electrical activity EAdi and the 
diaphragm pressure-time product PTP d j plotted over time during the volume 
maneuver (circles) and the two expulsive maneuvers at end-expiratory lung 
volume (EELV), one targeting a lower transdiaphragmatic pressure P di (squares) 
and the other a higher P d j transdiaphragmatic pressure P di (triangles). Figure 4 

30 shows that, for the representative subject, the center frequency CFdi declines 
more rapidly and to a greater extent during the volume maneuver (circles), which 
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required a high diaphragm's electrical activity EAdi for a similar diaphragm 
pressure-time product PTPdi, compared to the expulsive lower-pressure 
maneuver (squares). 

5 Table 2 



Individual CFdi values observed at end of each maneuver 



Subject 


Volume 
maneuver 


Expulsive maneuver 
Low pressure 


Expulsive 
maneuver 
High pressure 


1 


68.3 ± 8.7 


81.3 ±7.9 


61.6 ± 5.6 


2 


67.4 ±4.5 


86.7 ±8.0 


70.5 ± 13.2 


3 


82.0 ± 12.3 


94.3 ± 7.8 


71 .5 ±5.0 


4 


72.8 ± 5.8 


89.4 ±6.1 


73.4 ± 3.1 


5 


80.4 ± 9.3 


104.1 ±4.4 


83.6 ± 5.4 


Mean+SD 


74.2±6.8 


91.2±8.6 


72.1+7.9 


Values are means for each subject for eac 


h of the maneuvers performed. 



In order to produce a similar drop in center frequency CF di during the 
expulsive maneuver at EELV as was observed during the volume maneuver, 
10 more than a two-fold increase in the target transdiaphragmatic pressure P d i was 
required. This was associated with an increase in diaphragm's electrical activity 
EAdi from 25±8 % to 44±9% of maximum. As can be seen in Figure 4, the rate 
of decline of the center frequency CFdi was similar for the volume maneuver 
(circles) and the expulsive higher-pressure maneuvers (triangles). 

15 

Presented in Table 3 are the values of transdiaphragmatic pressure P d j, 
diaphragm's electrical activity EAdi and center frequency CF dj for the test-retest 
of the volume maneuver. During the retest, subjects successfully targeted a 
transdiaphragmatic pressure P d j that was similar to that generated during the 
20 initial volume maneuver (ICC=0.95). The diaphragm's electrical activity EAdi was 
also similar (ICC=0.93) as was the drop in center frequency CF di (ICC=0.98). 
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Table 3 

5 Test-retest of the volume maneuver 





CFdi 0 


CFdi 


EAdi 


Pdi 




(Hz) 


(Hz) 


(% max) 


(cm H 2 0) 


Sub- 
ject 


Vol 1 




Vol 2 


Vol 1 




Vol 2 


Vol 1 




Vol 2 


Vol 1 




Vol 2 


1 


90.4 




93.1 


68.3 




66.2 


46.5 




38.1 


47.8 




44.6 


2 


102.5 




108.2 


67.4 




65.7 


42.4 




44.4 


38.6 




38.3 


3 


94.7 




100.3 


82.0 




84.7 


49.8 




55.0 


34.1 




30.7 


4 


99.7 




99.5 


72.6 




72.1 


51.3 




50.6 


36.5 




35.8 


Mean 


96.8 




100.3 


72.8 




71.9 


66.7 




65.0 


25.3 




29.5 


(±SD) 


(5.4) 




(6.2) 


(6.7) 




(8.8) 


(10.7) 




(11.9) 


(9.4) 




(7.0) 


ICC 




0.94 






0.98 






0.93 






0.95 





EAdi, diaphragm electrical activity calculated as root-mean-square; CFdi 0 , baseline center 
frequency determined during resting breathing; CFdi plateau value of the center frequency at 
the end of the volume maneuver; P d j, transdiaphragmatic pressure; Vol 1, first volume 
maneuver performed; Vol 2, second volume maneuver performed; ICC, interclass correlation 



1 0 coefficient. 

DISCUSSION 

The study evaluated intermittent static contractions maintained at two 
15 different lung volumes, in order to examine the effect of altered neuromechanical 
coupling and increased diaphragm electrical activation, on diaphragm 
sarcolemma excitability, assessed by changes in center frequency CF d j. It was 
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found that, for a given targeted diaphragm pressure-time product PTP di , the drop 
in center frequency CF d j was greater when the diaphragm's electrical activity 
EAdi was increased by neuromechanical uncoupling, suggesting that the level of 
muscle activation influences the center frequency CF di . 

5 

Studies on the canine diaphragm have demonstrated that changes in 
center frequency CF d j are associated with changes in the mean action potential 
conduction velocity (APCV) [38], confirming previous mathematical models [29]. 
During muscle contractions, both center frequency CF d j and muscle fiber APCV 

10 depend to a smaller extent on the cable properties of the fiber [38, 39], and to a 
larger extent on the muscle membrane excitability [17, 18, 29, 39]. The 
excitability of the muscle fiber membrane is dependent on the trans-membrane 
gradient of potassium, and with increased muscle activation, efflux of potassium 
increases. In order to defend the extra-cellular potassium concentration and 

15 hence, the membrane potential, the cell depends on the re-uptake of potassium, 
e.g. via the ATP (Adenosine Triphosphate) dependent sodium/potassium pump 
[12], and washout via the blood circulation [25], i.e. diffusion of potassium from 
the extra-cellular space into the blood stream. Regardless if blood flow is 
reduced [23, 31, 42], or the muscle activation is increased, as in the present 

20 work, the muscles 1 electrical activity will indicate reduced membrane excitability, 
by shifts in the power spectrum toward lower frequencies. The center frequency 
CF di can also be affected by factors such as motor unit territory, number of fibers 
in the motor unit, dispersion in arrival times of the single contributions in the 
motor unit signal, dispersion in action potential conduction velocities between 

25 motor units that can cause the diaphragm's electrical activity EAdi power 
spectrum to shift [4, 29]. However, given that these influences are minor in 
healthy muscles [30] and given that the test situation did not allow for much 
variability in the contractile pattern, it is unlikely that these influences had more 
than a minor impact on the results. 

30 

In the study, a constant transdiaphragmatic pressure P di was targeted 
with a constant duty cycle at two different lung volumes, and it was therefore 
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assumed that transdiaphragmatic pressure P d i hindrance to blood flow under 
those conditions remained relatively similar at the different muscle lengths [19]. 
However, in order to achieve the same target transdiaphragmatic pressure Pdi at 
an increased lung volume, diaphragm's electrical activity EAdi was increased, 
5 which represents an increase in energy demand/consumption as well as 
increased metabolic output (e.g. potassium efflux) from the cell. As can be seen 
in Figure 3, the rate of decline of CF d i at increased lung volume was significantly 
higher than that observed when the same pressure was targeted at FRC 
(doubling of transdiaphragmatic pressure Pdi at the same lung volume) with 

10 lower diaphragm's electrical activity EAdi. Vitro studies have also demonstrated 
that increased activation (i.e. demand), accomplished by increasing stimulation 
frequency of a muscle shortened to 70% of its optimum length, in order to obtain 
the same tension generated at optimum length, resulted in an increased fatigue 
in the shortened muscle [14]. The current study therefore demonstrates that the 

15 higher diaphragm activation required for generating the target 
transdiaphragmatic pressure P di at an increased lung volume (neuromechanical 
uncoupling) influences the rate/extent to which center frequency CF d i decays. 
Further theoretical evidence for the impact of neuromechanical uncoupling on 
the center frequency CF d i is provided in the following description. 

20 

In the absence of neuromechanical uncoupling, an increase in 
transdiaphragmatic pressure P di is always associated with an increase in 
diaphragm's electrical activity EAdi. In the above reported study, doubling of 
transdiaphragmatic pressure P di at the same lung volume (FRC) was associated 

25 with an increase in diaphragm's electrical activity EAdi from 25% to 44% of the 
maximum. Beck et al [6] showed that diaphragm's electrical activity EAdi in 
absolute values is closely related to transdiaphragmatic pressure P dj , such that 
activation increases (i.e. energy demand increases) when pressure increases 
(i.e. energy supply decreases). However, this relationship is altered when the 

30 muscle length changes. In such a circumstance, the transdiaphragmatic 
pressure P d i continues to reflect diaphragm's electrical activity EAdi only when 
the transdiaphragmatic pressure P d j is normalized to the maximum 
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transdiaphragmatic pressure Pdi obtained at each corresponding lung volume 
[6]. It was previously shown that when the same diaphragm's electrical activity 
EAdi is targeted at different lung volumes, the higher resulting 
transdiaphragmatic pressure P d i generated at FRC promotes a greater drop in 
5 center frequency CF di than does the lower pressure produced at the higher lung 
volume [42]. Such results indicate that for a given neural activation, an increase 
in force or transdiaphragmatic pressure P d j reduces diaphragm excitability. 
Consequently, the use of the TT d j and pressure-time product as indices for 
predicting changes in the excitability of the diaphragm sarcolemma (as reflected 
10 by center frequency CF d j) is limited to conditions of constant neuromechanical 
coupling, where the diaphragm force generating capacity remains unaltered. 



Consistent with previous studies [3, 16, 21, 28], doubling of the target 
transdiaphragmatic pressure P dj at FRC in the present study increased the rate 
15 of decline of the center frequency CF di as well as the level to which it declined 
(Figures 3 and 4). This is partially explained by the increase in diaphragm's 
electrical activity EAdi, as discussed above. However, it is also partially 
explained by the fact that: 

20 i) diaphragm contractions with a higher transdiaphragmatic pressure P di 

tend to hinder blood flow (i.e. energy supply) relatively more than 
contractions producing a lower transdiaphragmatic pressure P di [19]; 
and 



25 ii) impaired blood flow to a muscle has the propensity to promote shifts 

in the electromyographic power spectrum toward lower frequencies 
[22, 30]. 



30 



Methodological and technical aspects 

In the study the contraction and relaxation periods were maintained at a 
fixed duration and therefore any potential influence of duty cycle on muscle 
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function [2, 22] was controlled for. It must be emphasized that accurate 
physiological measurement of the center frequency CF di depends on being able 
to control for: 

5 (a) changes in muscle-to-electrode distance; 

(b) electrode positioning with respect to the muscle fiber direction and 
location; 

10 (c) electrode configuration; 

(d) signal to noise ratio; 

(e) influence of cross-talk from other muscles (including the heart and the 
1 5 esophagus); and 

(f) electrode movement-induced artifacts [7, 36, 38, 39, 40]. 

In the study, the technology used to measure the power spectrum of the 
20 diaphragm's electrical activity EAdi spectrum included means for minimizing 
these influences [1, 36, 40]. The findings that evoked muscle action potentials 
are influenced by changes in lung volume [5, 15] have contributed to the 
assumption of a potential-inherent inaccuracy of measured amplitudes of the 
diaphragm's electrical activity EAdi [5, 15] and the center frequency CF d , [5]. 
25 However, during mild voluntary muscle contractions that do not alter diaphragm 
membrane excitability, it has been shown that chest wall configuration/lung 
volume and changes in muscle length have no effect on diaphragm's electrical 
activity EAdi and center frequency CF di [5, 6, 7, 17, 39]. Therefore the above- 
discussed effect of chest wall configuration/lung volume likely did not have an 
30 impact on the results. 



Another factor that could have influenced the results of the study is the 
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difference in partitioning the esophageal and gastric pressures for the same 
transdiaphragmatic pressure P d j during the various maneuvers. In a previous 
study [42], where subjects targeted the same diaphragm's electrical activity EAdi 
at higher and lower lung volumes, greater decreases in center frequency CF d j 
were consistently observed at EELV (higher transdiaphragmatic pressure Pen), 
regardless of whether subjects performed an expulsive (i.e. transdiaphragmatic 
pressure P di generated mainly by gastric pressure) or a Mueller maneuver (i.e. 
transdiaphragmatic pressure P di generated mainly by esophageal pressure) at 
EELV [42]. In a pilot trial to that study (unpublished observations), it was found 
that diaphragm contractions generating identical transdiaphragmatic pressure 
P dj , duty cycle and diaphragm's electrical activity EAdi, produced the same 
trajectory of decrease in center frequency CF di , whether subjects performed 
expulsive or Mueller maneuvers. Therefore, it is not believed that differences in 
the partitioning of the esophageal and gastric pressures during the volume and 
expulsive maneuvers in the current study had an effect on the outcomes 
observed. 

Clinical implications 

20 The results of the above reported study have direct implications to 

subjects or patients being weaned from mechanical ventilation. It is well known 
that patients undergoing a weaning trial may demonstrate dynamic changes in 
EELV (dynamic hyperinflation) [43], which similar to the study would alter the 
neuromechanical coupling of the diaphragm. In order to compensate for this 

25 uncoupling (i.e. maintain the same transdiaphragmatic pressure P d j), the patient 
would need to increase diaphragm activation. The combination of an increased 
activation of the patient's diaphragm, with an elevated transdiaphragmatic 
pressure P di would, according to the present study, lead to decreased center 
frequency CF di (excitability), and possibly an increased respiratory effort 

30 sensation [42]. Shifts in the H/L ratio of the power spectrum of the diaphragm's 
electrical activity EAdi have been reported in patients with respiratory failure in 
whom ventilatory assistance is removed [8, 13]. However, given that diaphragm 



10 
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weakness is prevalent in mechanically ventilated patients [24], it remains to be 
determined what combined levels of diaphragm's electrical activity EAdi and 
transdiaphragmatic pressure P dj would affect center frequency CF dj . 

5 % CONCLUSION 

The above-reported study shows that diaphragm activation can be used 
to determine diaphragm membrane excitability and changes in center frequency 
CF d j. Furthermore it shows that the diaphragm pressure-time product PTP dl and 
10 tension-time index TT di cannot be considered as valid reflections of diaphragm 
energy consumption and/or sarcolemma excitability when neuromechanical 
coupling is altered. 

With data from the above investigation or study, the diaphragmatic 
1 5 muscle force can be estimated from measurements of the diaphragm's electrical 
activity EAdi in two ways. 

A first way for estimating the diaphragmatic muscle force uses the 
following equation: 



20 



F = juEAdi (1) 



where F is the diaphragmatic muscle force, |n is a proportionality constant, and 
EAdi is a measure of the signal strength of the electrical activity of the patient's 
25 diaphragm. Here the square root of the first power spectral moment is used 
since it represents the signal strength, which has been compensated for the 
influence of changes in the propagation velocity of the myoelectric action 
potentials [29]. 



30 



A second way for estimating the diaphragmatic muscle force uses the 
spectral changes during diaphragm contraction. For a forceful periodic muscle 
loading, the center frequency CFdi decreases from an initial center frequency 
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CF 0 to a final plateau value CFoo according to the equation [26]: 



CF„ = CF Q (1 - k)T d /[(l - k)T d + kT r ] (2) 

5 where k is the duty cycle, i.e. the inspiration time in relation to the total time 
period, and T R is the center frequency CF d j recovery time constant pertaining to 
an approximately exponential time curve which is rather independent of the 
muscle force [9]. The symbol T D denotes the time constant for the decrease in 
center frequency CF di , which is related to the muscle force as [27]: 

10 

T D ^?ll{F-F c ) (3) 



In this equation ti is a proportionality constant and F c is a critical force level 
above which muscle fatigue starts to develop. Equation (2) is rearranged to 
1 5 obtain the experimentally determinable quantity: 

Q = T r /T d =[(1-k)/k] [(CF Q -CFJ/CFj (4) 



20 



Equations (3) and (4) then give: 

F = F c +Qr,IT R (5) 



Making equal the two force estimates of equations (1 ) and (5) the 
following relation is found: 

25 

aEAdi-j3~Q = 0 (6) 



where 



30 



a = juT R /?j 
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and 

P = F c T R ir, (8) 

5 

Relation (6) represents a set of three equations (for the three experimental 
conditions) with two unknowns. A fitting procedure with data from the following 
Table 4 with simultaneous minimization of the relative errors in the diaphragm's 
electrical activity signal strength EAdi and the quantity Q, gives the values a = 
1 0 0.0041 7 and p = 0.041 9 with a relative fitting error of 0.24. 



Table 4 

Experimental results and calculated values 















Force ratios 


Geometric factors 




Pdi 
(cm H 2 0) 


EAdi 
(a.u.) 


CF 0 
(Hz) 


CFoo 
(Hz) 


Q 




<l>ii 




Yi 

(cm H 2 0) 


Yii 

(cm H 2 0) 




Volume Maneuver 


31.2 


60.9 


100 


74.2 


0.175 


6.12 


5.20 


5.66 


122 


144 


133 


Expulsive maneuver 
(Lower pressure) 


31.6 


24.9 


100 


91.2 


0.050 


2.50 


2.19 


2.35 


303 


347 


325 


Expulsive maneuver 
(Higher pressure) 


67.5 


44.3 


100 


72.1 


0.196 


4.45 


5.69 


5.07 


364 


285 


324 



P di , transdiaphragmatic pressure; a.u., arbitrary units; EAdi, signal strength of the electrical 



1 5 activity of the diaphragm; CFdi 0 , baseline diaphragm center frequency determined during resting 
breathing; CFdi, plateau value of the diaphragm center frequency at the end of the maneuver; Q, 
ratio of the time constants of CFdi recovery and decline, see equation (5); <)>,, see equation (9); <j>n, 
see equation (10); ty m , mean of <|>| and 4>n; Yi« see equation (14); y iU see equation (15);y mi mean of yi 
and Yij. 

20 



With a and p known, the experimental values of the diaphragmatic 
muscle force F can be expressed in relation to the critical force level Fc for onset 
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of deterioration of cell excitability, i.e. the critical force level above which muscle 
fatigue starts to develop. The two ways to describe this are obtained by 
rearranging equations (1) and (7), and equation (5), respectively, which gives: 

<f>j =(F/F c ) f = aEAdi/j3 (9) 



and 



10 



30 



* n =iFIF c ) B =\ + Qifi 



(10) 



These quantities have been determined and are listed in Table 4 together with 
their mean values (j) m . 

The observed transdiaphragmatic pressure P dl is assumed to be related 
15 to the diaphragmatic muscle force F as: 

Pd/ = FG (11) 

where G is a geometrical factor taking into account that the diaphragm muscle 
20 changes its shape with the inspired volume. This factor G is thus assumed to be 
the same during the expulsive maneuvers with lower or higher P d j production 
performed at end-expiratory lung volume. As with the force relations, the 
transdiaphragmatic pressure P d j can be expressed in two ways, relating to the 
diaphragm's electrical activity signal strength EAdi and to the fatigue induced 
25 spectral changes. Combining equations (1), (5), and (11) leads to the following 
relations: 

Pdi = ju EAdi G (12) 

and 



Pdi = (F c +Q? 7 /T R )G (13) 
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Relations (12) and (13) can be further developed with relations (7) and 
(8) into the two following relations: 

5 Yi =(GT}/T R )j=Pdi/(a EAdi) (14) 

and 

(15) 

10 Numerical values, calculated for the two expressions, are given in Table 4 
together with their mean values y m . 

From the results listed in Table 4, it can be concluded that the 
diaphragmatic muscle force F in relation to the critical force level F c are 

15 approximately the same during the volume maneuver and the higher pressure 
expulsive maneuver, which is also reflected in their deterioration of cell 
excitability, expressed by the factor Q. During all conditions the diaphragmatic 
muscle forces F are above the critical force level F c as shown by values of § m in 
Table 4. The geometrical dependence, expressed by the factor ym, is obviously 

20 the same during lower pressure expulsive maneuver and higher pressure 
expulsive maneuver, but is much less during the volume maneuver. The ratio 
between the y values in the volume maneuver and the expulsive maneuvers is 
about 0.41. Since the r| values and the T R values are expected to be 
independent of the maneuvers, this means that also the factors G have the 

25 same ratio. This indicates a much lower efficiency to convert force into pressure 
during the volume maneuver. The tension time index TT dj , taking into account 
the timing and the pressure, is thus not sufficient to describe the complexity of 
the fatigue development. At least it has to be modified with a volume dependent 
correction factor. Better, though, are methods reflecting the deterioration of cell 

30 excitability and not the mechanical result of the contraction. 
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Electromyographic and mechanical methods to detect muscle fatigue 

Based on the above results, techniques to determine critical levels of 
muscle fatigue during periodic loading (such as respiration) will be described. A 
5 number of equations relating certain physiological variables to each other are 
needed and they will be derived prior to the description of these techniques. 

Periodic muscle load characteristics 

10 Consider a periodic muscle loading, such as the respiratory work, in 

which repeated muscle contractions alternate with muscle relaxations. The 
periodic muscle loading is characterized by a time period T 0 and its two parts: 
the duration of muscle contraction Ti and the duration of muscle relaxation T 2 
where: 

15 

To = Tn + T 2 (16) 

In order to simplify the equations, the duty cycle k is determined as: 

20 k = T1IT0 (17) 

The mean diaphragmatic muscle force developed during the time interval 
T1 is denoted F. 

25 Myoelectric changes due to fatigue 

Isometric fatiguing contractions cause the center frequency CFdi of the 
diaphragm's electrical activity EAdi diaphragm's electrical activity to decrease 
exponentially from its resting value CF 0 with a time constant Tf- During recovery 
30 the center frequency CFdi returns gradually to its normal value following an 
approximately exponential course, described by the recovery time constant T R . It 
is observed that many other characteristics of the power spectrum of the 
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diaphragm's electrical activity EAdi exhibit the same dependencies such as the 
median frequency, the zero crossing density, the so-called hi-over-low value, 
etc. The recovery time constant depends mostly on the density of capillaries in 
the muscle and is rather insensitive to the exerted force. The fatigue time 
5 constant is strongly dependent on the force when it exceeds a certain critical 
level F c - The relation is: 

T F = /?/(F-Fc) for F>F C (18a) 

10 and 

T F -> oo for F<F C (18b) 

The combination of repeated work and recovery events causes the center 
15 frequency CFdi to decrease from the initial value to a final plateau value CFoo, at 
which there is a balance between the metabolite production during work and 
wash-out during recovery. The plateau value is: 

CFco = CF 0 (1 - at)7 f /[(1 -*)7> + kT r ] (19) 

20 

Introducing the notations: 

ACF = CFq-CFoo (20) 

25 and 

s = ACF/CFo (21) 



30 



Equation (19) can then be rearranged to read: 
k = 1 /[1 +(T R /T F )(ACF/ OF,,)] 



WO 2005/077268 



PCT/CA2005/000217 



With the notation: 



Q = T R IT F 



(23) 



5 



it is found that: 



Q = [(1 - k) I k] ACF I CF C 



CO 



(24) 



which is an experimentally measurable quantity. 



10 



Force and pressure 

The diaphragmatic muscle force F can be determined for skeletal 
muscles working over joints without synergistic effects from other muscles. For 
15 the diaphragm muscle the force cannot be directly measured, rather the 
transdiaphragmatic pressure P d i is obtained as a proportional measure. The 
following relation could be used: 



where \x is a proportionality constant and E is the signal strength of the 
diaphragm's electrical activity EAdi, preferably based on the first spectral 
moment which is rather insensitive to metabolic changes caused by fatigue. The 
relation to the pressure is proportional but non-linear. This fact is taken into 
25 consideration by introducing the factor G(V) which is volume (V) dependent, i.e.: 



F = juE 



(25) 



20 



Pdi 



= FG(V) 



(26) 



Thus, 



30 



juG(V) = P di /E 



(27) 
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which also is an experimentally measurable quantity. 

Myoelectric signal strength and spectral changes 

5 Rearrangement of equation (18a) and insertion of equations (23) and (25) 

gives: 

ccE - J3 - Q = 0 (28) 

10 where 

a = juTr/tj (29) 

and 

15 

P = FoTr/tj (30) 

It can be observed that a is dependent, through the parameter jn, on the 
electrode geometry and placement in relation to the muscle, while the other 
20 parameters are rather constant for similar muscles. 

Experiments under fatiguing conditions at any volume give corresponding 
values of E and Q (through the center frequency changes). A data fitting 
procedure (not regression) gives numerical values to a and p. With a and p 
25 known, an estimate of the diaphragmatic muscle force F can be obtained in 
relation to its fatigue threshold value, i.e.: 

F/F c = Ea/jS (31) 

30 As long as F/F c is smaller than one, isometric fatigue of the patient's 

muscle does not develop. That means that the signal strength should be lower 
than the critical value: 
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E < £,som = P/a (32) 

\ 

For periodic muscle work, higher forces and signal levels are tolerable. 

5 

Spectral changes as indicators of tolerable concentration of metabolites 

The relative spectral change 8 of the diaphragm's electrical activity EAdi, 
defined in equation (21), is an indirect measurement of remaining concentration 

10 of metabolites in the muscle during periodic fatiguing contractions. It seems that 
the muscle very rapidly goes into an anaerobic metabolic state once the force is 
higher than F c and that virtually all contractions above this level causes changes 
in the the center frequency CFdi. Therefore it is likely that a certain small value 
of 8 is tolerable as long as it is below a certain critical level, which we denote e c . 

15 With this critical value introduced into equation (22) and simultaneous use of 
equations (18a) and (23), it can be found that a condition for long term fatigue 
not to occur is: 

k< 1/{1 +[(1 -s c )l e c ]T R {F-F c )l r!} (33) 

20 

This expression can be rearranged to give the force condition: 

F < F c +[(1-*:)/ir][fl b /(1-*)]j7/r /? (34) 
25 or, together with equation (30), 

F < F c {1 + [(1 -k)/k][8c/(1 -ec)]/P> (35) 

Since the force in diaphragmatic contractions cannot be simply measured, 
30 equations (33) to (35) are expressed as functions of the signal strength E and 
the transdiaphragmatic pressure Pdi. Use of equations (25) and (26) give for the 
signal strength E of the diaphragm's electrical activity EAdi: 
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5 



k < 1/{1 +[(1 -ec)/ec](aE-p)} (36) 

and 

E < {p + [(1 -k)/k][8c/(1 -e c )]}/a (37) 
and for the transdiaphragmatic pressure Pdi: 
10 k < 1/{1 +[(1-ec)/ec](aP di -p)} (38) 

and 

Pdi < (xG(V){p + [(1 -k)/k][sc/(1 -ec)]} (39) 

15 

NON-RESTRICTIVE ILLUSTRATIVE EMBODIMENT OF A METHOD 
AND DEVICE FOR DETERMINING AN OPTIMAL LEVEL OF 
VENTILATORY ASSIST TO A VENTILATOR DEPENDENT PATIENT 

20 

Non-restrictive illustrative embodiments of the method and device for 
determining an optimal level of ventilatory assist to a ventilator-dependent 
patient will now be described. 

25 First embodiment of Figure 5 

Operation 501 

The signal strength of the diaphragm's electrical activity EAdi is 
30 monitored through a detector 502. As illustrated in Figure 1, detector 502 may 
comprise, for example, a computer 4 to measure the signal strength of the 
electrical activity EAdi of the patient's diaphragm through a linear array 5 of 
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electrodes mounted on an esophageal catheter 6 inserted through the patient's 
nostril (or patient's mouth) until the electrode array 5 is positioned in the gastro- 
esophageal junction 10 of the patient's diaphragm 7. 

5 

Operation 502 

A calculator 503 calculates the coefficients a and p using equation (28): 
10 aE - p - Q = 0 (28) 

with myoelectric data from fatigue tests (calibration). Fatigue test can be 
performed by either reducing the level of assist, or performing a short airway 
occlusion, while measuring the myoelectric activity during a few inspiratory 
15 attempts. To shorten and facilitate the fatigue test the subject could be 
encouraged to voluntarily increase his efforts. Such a test is routinely performed 
to determine the maximum inspiratory airway pressure. 

Operation 504 

20 

The calculator 503 calculates the duty cycle k as described hereinabove. 
Operation 505 

25 The calculator 503 calculates estimates of a critical level of the relative 

spectral change e c of the diaphragm's electrical activity EAdi from the general 
experimental fact that fatigue does not occur below a duty cycle of 0.2 even at 
maximum muscle force and that the critical force level F c is approximately 0.2 
times the maximum force. Equation (33) then gives s c & p/((3+8/9), or, since both 

30 8 C and p are small quantities: 



e c & J3 



(40) 
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Operation 506 

The calculator 503 calculates a critical signal strength of the diaphragm's 
5 electrical activity EAdi above which isometric muscle fatigue develops, using the 
relatiion: 

E < Eisom = P/oc (32) 

10 

Operation 507 

If myoelectric monitoring is used (giving signal strength and duty cycle), 
the calculator 503 calculates a critical signal strength of the diaphragm's 
15 electrical activity EAdi above which long term muscle fatigue develops, is 
calculated using equation (37): 

E < {p + [(1-K)/K][ec/(1-ec)]}/a (37) 

20 Operation 508 

A controller 509 controls the ventilatory assist, for example the gain of the 
ventilatory assist at a level such that the signal strength of the diaphragm's 
electrical activity EAdi does not exceed that described in relation (37) (higher 
25 support suggest unnecessary muscle inactivation) to prevent long-term muscle 
fatigue to develop:: 

E < {(3 + [(1-K)/K][ec/(1-ec)]}/a (37) 

30 However, the signal strength of the diaphragm's electrical activity EAdi 

should not exceed that described in equation (32) (this level indicates the level 
for muscle fatigue during isometric contractions) to prevent isometric muscle 
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fatigue to develop:: 

E < E ISO m = j3/a (32) 
5 Second embodiment of Figure 6 
Operation 601 

The signal strength of the diaphragm's electrical activity EAdi is 
10 monitored through a detector 602. As illustrated in Figure 1, detector 602 may 
comprise, for example, a computer 4 to measure the signal strength of the 
electrical activity EAdi of the patient's diaphragm through a linear array 5 of 
electrodes mounted on an esophageal catheter 6 inserted through the patient's 
nostril (or patient's mouth) until the electrode array 5 is positioned in the gastro- 
1 5 esophageal junction 1 0 of the patient's diaphragm 7. 

Operation 603 

A detector 604 monitors the patient's transdiaphragmatic pressure P dj . As 
20 illustrated in Figure 1, detector 604 may comprise, for example, a computer 4 to 
continuously measure the transdiaphragmatic pressure P d i by detecting the 
esophageal P es and gastric P ga pressures through respective gastric 8 and 
esophageal 9 balloons mounted on the catheter 6 on opposite sides of the array 
5 of electrodes, and by processing the detected esophageal P es and gastric P ga 
25 pressures to obtain the patient's transdiaphragmatic pressure P d j. 

Operation 605 



30 



A calculator 606 calculates the coefficients a and p using equation (28): 
<xE - J3 - Q = 0 (28) 
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with myoelectric data from fatigue tests (calibration). Fatigue test can be 
performed by either reducing the level of assist, or performing a short airway 
occlusion, while measuring the myoelectric activity during a few inspiratory 
attempts. To shorten and facilitate the fatigue test the subject could be 
5 encouraged to voluntarily increase his efforts. Such a test is routinely performed 
to determine the maximum inspiratory airway pressure. 

Operation 607 

1 0 The calculator 606 calculates the duty cycle k as described hereinabove. 

Operation 608 

The calculator 606 calculates estimates of a critical level of the relative 
15 spectral change e c of the diaphragm's electrical activity EAdi from the general 
experimental fact that fatigue does not occur below a duty cycle of 0.2 even at 
maximum muscle force and that the critical force level F c is approximately 0.2 
times the maximum force. Equation (33) then gives s c ~ f3/(j3+8/9), or, since both 
£ C and p are small quantities: 

20 

sc * P (40) 
Operation 609 

25 The calculator 606 calculates a critical signal strength of the diaphragm's 

electrical activity EAdi above which isometric muscle fatigue develops, using the 
relation: 



30 



E < Eisom = fi/a 
Operation 610 



(32) 
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The calculator 606 calculates a critical level of the transdiaphragmatic 
pressure P d j using relation (39): 

P di < jiG(V){p + [(1-K)/K][sc/(1-ec)]} (39) 

5 

Operation 610 requires knowledge about the geometrical G(V) 
dependence. This factor G(V) can be obtained from a calibration of the 
experimentally measurable quantity P d j / E as shown in equation (27): 

10 juG(V) = P di /E (27) 

Alternatively, geometrical dependence G(V) of inspiratory pressure can 
also be estimated by performing single or multiple breath airway occlusions at 
two lung volumes, e.g. end-inspiration and end expiration lung volumes, while 
15 the volume difference is measured by a computer with, for example, at least one 
flow meter (see computer 4 and pneumotachograph 3 of Figure 1 ). 

Operation 611 

20 A controller 61 2 controls the ventilatory assist, for example the gain of the 

ventilatory assist at a level such that: 

- the monitored signal strength of the diaphragm's electrical activity EAdi does 
not exceed that described in relation (32) (this level indicates the level for 
25 muscle fatigue during isometric contractions) to prevent isometric muscle fatigue 
to develop: 

E < Eisom = P/oc (32) 



30 



- the monitored patient's transdiaphragmatic pressure Pdi does not exceed that 
described in relation (39) to prevent long-term muscle fatigue to develop: 
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Pdi < MG(V){p + [(1 -k)/k][6c/(1 -8c)]} (39) 

Although the present invention has been described hereinabove with 
reference to non-restrictive illustrative embodiments thereof, it should be kept in 
5 mind that these embodiments can be modified at will within the scope of the 
appended claims without departing from the spirit and nature of the present 
invention. In particular but not exclusively: 

- the present invention pertains not only to CFdi and RMS but possibly to other 
1 0 types of measures; 

- the present invention can be implemented through measurement of the 
electrical activity of respiration-related muscles other than the diaphragm; 
and 

15 

- the present invention is concerned with any method of mechanical 
ventilation, including negative pressure ventilation. 
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